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Abstract 

We investigate dynamically and statistically diffusive motion in a chain of linearly coupled 
2-dimensional symplectic McMillan maps and find evidence of subdiffusion in weakly and 
strongly chaotic regimes when all maps of the chain possess a saddle point at the origin and 
the central map is initially excited. In the case of weak coupling, there is either absence of 
diffusion or subdiffusion with q > 1-Gaussian probability distributions, characterizing weak 
chaos. However, for large enough coupling and already moderate number of maps, the system 
exhibits strongly chaotic [q « 1) subdiffusive behavior, reminiscent of the subdiffusive energy 
spreading observed in a disordered Klein-Gordon Hamiltonian. Our results provide evidence 
that coupled symplectic maps can exhibit physical properties similar to those of disordered 
Hamiltonian systems, even though the local dynamics in the two cases is significantly different. 

Keywords: Complex statistics, multi-dimensional maps, McMillan map, Klein-Gordon disordered 
Hamiltonian, chaotic and diffusive motion, q-Gaussians, Tsallis entropy 

In memory of our beloved friend and distinguished scientist Professor Theodore Papatheodorou, 
who taught us so eloquently that numerical analysis is indispensable when one tries to penetrate the 
mysteries of realistic physical systems. 


1 Introduction 

The absence of diffusion in disordered media (the so-called Anderson localization [2]) is an impor¬ 
tant effect that arises in transport phenomena, such as electromagnetic, acoustic and spin waves, 
in different types of classical or quantum linear systems. It is more interesting, however, to study 
what happens to the disordered system if nonlinearity is introduced. Understanding the effect of 
nonlinearity on the localization properties of wavepackets in disordered systems has attracted the 

attention of many researchers to date [HI [TH) (TU [Ml [371 [311 ESI (HI 13111131 Ull Hi HI 113] ■ Most 

of these works consider the evolution of an initially localized wavepacket in a chain of particles 
modeled by a continuous dynamical system and show that it spreads subdiffusively in moderately 
nonlinear systems, while for strong enough nonlinearities a substantial part of the wavepacket re¬ 
mains self-trapped. In particular, for single-site excitations, the spreading of the wavepacket leads 
to an increase of the second moment m 2 in time of the form m 2 ^ both in the diffusive and 
self-trapping regimes [Hllia 1331137]. 

Currently, a greatly debatable question concerns the long time behavior of wavepacket spread¬ 
ing in disordered nonlinear particle lattices. Recently, it was conjectured [BIS] that chaotically 
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spreading wavepackets will asymptotically approach Kolmogorov-Arnold-Moser (KAM) torus-like 
structures in phase-space, while in all numerical simulations chaotic spreading shows no sign of 
slowing down [111111133]. Nevertheless, for particular disordered nonlinear models, some numerical 
indications of a possible slowing down of spreading have been reported |25l [24] . In Ref. [3] , employ¬ 
ing ideas from nonextensive statistical mechanics |35] , we studied the same problem and found no 
such sign of quasi-periodic KAM behavior and all our evidence appears to strengthen the conjecture 
that waves spread subdiffusively and chaotically for arbitrarily long times in nonlinear disordered 
media. It would be interesting, therefore, if similar results could also be found in coupled nonlinear 
maps modeling nonlinear disordered media as they are computationally easier to study than the 
systems of ordinary differential equations of their Hamiltonian counterparts. Hamiltonian systems 
and symplectic maps, and in general discrete time systems, can be identified through discretization 
techniques, i.e. Poincare sections. The drawback of this approach however is that the resulting 
map is in general, not explicitly given and thus one has to resort to the construction of a suitable 
discrete system that is able to exhibit similar dynamical and statistical properties as the disordered 
Hamiltonian counterpart. 

In the present paper, we use tools of nonextensive statistical mechanics to investigate the 
connection between regimes of “weak” and “strong” chaos and diffusion in a system of N coupled 
symplectic McMillan maps (CMM) [331113], with the purpose of linking them to analogous results 
obtained recently for a Klein-Gordon (KG) Hamiltonian system in the presence of disorder and 
nonlinearity [3]. In this framework, we are interested to explore the question whether coupled 
symplectic maps can be used to study subdiffusive spreading of wavepackets in disordered media, 
instead of using continuous flow dynamical models such as the KG Hamiltonian. 

As is well-known, probability distribution functions (pdfs) of chaotic trajectories of dynamical 
systems have been studied by many authors, aiming to understand the transition from deterministic 
to stochastic dynamics [31|29l[l2]. One of the most fundamental questions concerns the existence 
of an appropriate invariant probability density or ergodic measure, characterizing chaotic motion 
in phase-space. If it is possible to demonstrate the existence of such a measure, one can postulate 
that the system is at thermal equilibrium and can therefore be studied from the point of view of 
classical statistical mechanics. 

Since the existence of an invariant measure is not known a priori, one can still proceed in the 
context of the central limit theorem m and consider the variables of a chaotic solution at discrete 
times ti, z = I,..., M as realizations of Af independent and identically distributed (iid) random 
variables j = 1,... ,Af. If the motion is uniformly ergodic in some region of phase-space, 

one finds that the pdfs of the sums of these variables converge rapidly to a Gaussian distribution, 
whose mean and variance are those of the In such cases we call the dynamics “strongly” 

chaotic, since at least one Lyapunov exponent is positive [61IZ1I3Q] and the respective subset of the 
constant energy manifold is uniformly covered by chaotic orbits, for all but a (Lebesgue) measure 
zero set of initial conditions. 

Now, if the system has not yet reached thermal equilibrium, it may happen that its orbits “stick” 
for long times on the boundaries of islands surrounding stable periodic orbits, where Lyapunov 
exponents [HE] become very small and may even vanish [301135] . In such regimes, the motion 
may be called “weakly” chaotic, as trajectories get trapped within complicated sets of cantori and 
diffuse through multiply connected domains of phase-space in a highly nonuniform way 
nnEn]. Many such examples have been discovered in physically realistic systems studied recently 

[311 [HI [13133] ■ 

In the present paper, we compare the diffusive dynamics of a system of coupled symplectic 
McMillan maps [20l [15] , with analogous results obtained recently in Ref. [3] for a KG Hamiltonian 
1-dimensional lattice in the presence of disorder and nonlinearity. In particular, we find in both 
models that weakly chaotic dynamics is associated with subdiffusive motion and is characterized 
by pdfs of sums of variables that do not rapidly converge to a Gaussian, but are well approximated 
for long times by the so-called g-Gaussian distribution [33] 


P(s) 


aexp (—/3s^) = a 


1 - (1 - q)l3s" 


( 1 ) 


where the q entropic index satisfies 1 < g < 3, /3 is an arbitrary parameter and a is a normalization 
constant. Eventually, of course, it is expected that chaotic orbits will wander within larger chaotic 
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seas, where the dynamics is more strongly ergodic. When this occurs, the q entropic index of 
Eq. O decreases towards 9 = 1, which represents the limit at which the pdf becomes a Gaussian 
distribution (see Fig. [3] for an illustrative example). 

In Ref. [3], we studied wavepacket spreading in a Hamiltonian system represented by a dis¬ 
ordered nonlinear KG chain of = 1000 particles liiini and found that even though there are 
phase-space regions of weak chaos, strongly chaotic dynamics eventually prevails characterized by 
q-Gaussian pdfs that approach a Gaussian pdf for long times. We have thus concluded that the 
motion will never approach a KAM regime of invariant tori as suggested by some authors [HE]. 

In analogy with these findings, the statistics of diffusive motion in chains of N CMM maps 
yield similar results to those for the KG system, although there are substantial differences in the 
dynamics of the two systems. In particular, in the KG system the dynamics is such that each 
particle has an elliptic point at the origin, while in our study for the GMM maps each map has a 
saddle point at the origin. Moreover, the KG system has a conserved quantity (see Eq. ([9|)), while 
there is no such known quantity for the studied CMM maps of Eqs. m- Another difference is that 
the KG system is characterized by diffusive processes of one epoch m whereas in our study we hnd 
diffusion described by two different epochs. Finally, in the KG system large nonlinearities result in 
selftrapping, while smaller nonlinearities are connected to subdiffusion. In contrast, in the CMM 
maps large coupling strengths lead to diffusion while smaller values show no spreading. Despite 
these differences, we do find evidence of subdiffusion accompanied by weak and strong chaos, as 
in the case of the disordered KG chain, when the maps have a saddle point at their origin and the 
coupling strength is high enough. In that sense our results demonstrate that coupled symplectic 
maps can be used to model diffusive motion. 

Our paper is organized as follows: In Sec. [2] we outline the details of our study of statistical 
distributions corresponding to weakly and strongly chaotic behavior, while in Sec. |3]we summarize 
the results of Ref. [3] to help the reader compare with our results on the CMM system presented 
in Sec. m Finally, Sec. El contains our conclusions and discussion. 


2 Computation of statistical distributions of weak and strong 
chaos 

In this section we focus our attention on the statistical properties of chaotic diffusion in a 2N- 
dimensional CMM system of N linearly coupled 2-dimensional symplectic maps evolving in discrete 
time n. We wish to explore diffusion in regimes of weak chaos, where Lyapunov exponents are 
positive but very small. Such situations often arise when solutions move slowly through thin 
chaotic layers, wandering through a complicated network of higher order resonances, often sticking 
for very long times to the boundaries of islands constituting the so-called “edge of chaos” regime 

[33]. 

Many interesting questions can be asked in this context: How long do these weakly chaotic 
states last? What types of pdfs characterize their statistical behavior and how can one relate them 
to the diffusive properties of the motion? Is there a connection between complex statistics and 
diffusion as the dimensionality of the system increases? What is the effect of different parameters 
and what role does disorder play? 

To answer these questions, we use the orbits of a multi-dimensional GMM model to construct 
pdfs of suitably rescaled sums of M values of a generic observable rji = i = 1,..., M that 

depends linearly on the coordinates x of the orbit. Viewing these as iid random variables in the 
limit of M —>■ 00 , we evaluate their sum 
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for j = 1,..., different initial conditions and study the statistics of Eq. 
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their mean value (S^j) = rescaled by their standard deviation ctm 
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Plotting the normalized histogram of the probabilities P(sj^^) as a function of we compare our 
pdfs with a g-Gaussian of the form 


P(sm ) = aexp ) = 


1 - (1 - q)l3s 


0)2 
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cf. ([!]), where q is the so-called entropic index. Note that this is a generalization of the well-known 
Gaussian pdf, since in the limit q —>■ 1 we have exp^{—f3x^) —>■ exp(—/3ai^). Moreover, it can be 
shown that the g-Gaussian distribution o is normalized when 
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where P is the Euler P function. Glearly, Eq. m shows that the allowed values of q are 1 < g < 3 
for this normalization. 

The index q appearing in Eq. o is connected with the Tsallis entropy |35| 


Sq = k 


9-1 


with 


w 

2=1 


(7) 


where i = 1,... ,W counts the microstates of the system, each occurring with a probability Vi and 
k is the well-known Boltzmann constant. Just as the Gaussian distribution represents an extremal 
of the Boltzmann-Gibbs entropy S'gQ = Si = kY^i Vi \n.Vi, so is the g-Gaussian ((T|) derived by 
optimizing the Tsallis entropy of Eq. 0 under appropriate constraints. 

Systems characterized by the Tsallis entropy are said to lie at the “edge of chaos” and are 
significantly different from Boltzmann-Gibbs systems, in the sense that their entropy is nonadditive 
and generally nonextensive [35]. In fact, a g-central limit theorem has been proved |36j for q- 
Gaussian distributions © that is of the same form as the classical central limit theorem. 

As we demonstrate in the following sections, in regions of weak chaos these distributions are 
well-fitted by a g-Gaussian pdf for fairly long time intervals, whose q value is distinctly greater 
than f. It may happen, of course, for longer times that the orbits begin to diffuse through domains 
of strong chaos, in which case q tends to 1 and the well-known form of a Gaussian pdf is eventually 
recovered. 

To study diffusive processes of chaotic dynamics in multi-dimensional CMMs, we adopt the 
classical definition of the mean square displacement (MSD) (r^) given by the normalized sum of 
the squares of the distances of the orbits from their starting point. If the motion is a random 
walk, we expect that 

(r^) = Dn\ (8) 

where n denotes the iterations of the dynamics, 7 < 1 characterizes the process as subdiffusive, 
7=1 corresponds to normal diffusion and 7 > 1 to super-diffusive motion. D is the diffusion 
coefficient and 7 = 2 corresponds to ballistic motion. Finally, we compute the corresponding 
Lyapunov exponents of the dynamics following the well-known methodology presented in Refs. 

laiTiisni- 
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3 Chaotic diffusion in a KG particle chain 

In the case of a high-dimensional KG disordered particle chain, one may concentrate on a low 
energy (subdiffusive) and a higher energy (self-trapping) case and verify that subdiffusive spreading 
always occurs following specific power-laws. Integrating the equations of motion for long times and 
computing probability distributions of sums of the positions of particles, we have found in Ref. 
[3] convincing evidence that the dynamics does not relax onto a quasi-periodic KAM torus-like 
structure, but continues to spread chaotically along the KG chain of particles for arbitrarily long 
times. 

To show this, we considered a quartic KG lattice described by the N degree of freedom Hamil¬ 
tonian 

Hkg = ^ ^ -I- — xiY = E, (9) 

1=1 

where xi and pi are the generalized positions and momenta on site I, and ii are chosen uniformly 
randomly from the interval [^, |] to account for the disorder present at each site I of the chain. 
This Hamiltonian conserves the value of the total energy A > 0 of the system, which, for fixed 
disorder strength W, serves as a control parameter of the nonlinearity. In our study, we followed 
the evolution of single site excitations solving the equations of motion 

3 1 

xt = -eixi - Xi + -^{xi+i + xi-i - 2xi), l = l,...,N, (10) 

with fixed boundary conditions (i.e. xq = xn+i = 0 ) and monitored the corresponding normalized 
energy density distributions. 

In Fig. [T] (taken as is from Ref. i) we exhibit two representative examples of numerical pdfs 
with different q entropic indices for a low energy subdiffusive case with E = 0.4 and N = 1000 (the 
initial value of a; 5 oo( 0 ) is adjusted so that E = 0.4). In panel a) we plot the numerical distribution 
(dashed curve) for the observable rji = X 500 computed in the time interval [0,10®] for M = 10^ and 
note that it is well fitted by a g-Gaussian distribution (solid thick curve) with qi = 0.993 ± 0.009. 
This is a case where the numerical distribution is indistinguishable from a Gaussian (q = 1) plotted 
as a dotted curve. On the other hand, panel b) which is the same plot as a) for the observable 
7729 = Sf=486 reveals a clear g-Gaussian distribution ([T|) with (729 = 1.22 ±0.01, over nearly four 
decades on the vertical axis. 




Figure 1: Characteristic pdfs for the disordered N — 1000 KG system for single site 
excitations. Panel a): Plot of the numerically computed pdf (dashed curve) for the observable 
Vi = 2:500 in the time interval [0,10®] with qi = 0.993 ±0.009, taken from fitting with a g-Gaussian 
distribution o in solid thick. Panel b): Similar plot of the numerically computed pdf (dashed) 
for the observable 7729 = ^ interval [ 0 , 10 ®] with <729 = 1-22 ± 0 . 01 , taken from 

fitting with a g-Gaussian distribution o in solid thick. In both panels, N = 1000, M = 10^ and 
E = 0.4 that corresponds to the subdiffusive case. Note that the vertical axes are in logarithmic 
scale, while the dotted curve is the Gaussian pdf (i.e. <7 = 1). 

Reproduced from Chaos 24, 024405 (2014). 

© 2014 AIP 
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In this way, we studied in Ref. the dynamics and statistics of diffusive motion in the 
1-dimensional nonlinear disordered KG particle chain and obtained results that can be directly 
compared to what we describe in the present paper using CMM maps. In particular, for a dis¬ 
ordered KG chain of TV = 1000 particles, we focused on a low energy (subdiffusive) and a higher 
energy (self-trapping) case and verified that subdiffusive spreading always occurs following specific 
power-laws with exponents smaller than 1 as pointed out in the literature M- Moreover, integrat¬ 
ing the equations of motion for times as long as 10® and computing the corresponding pdfs, we 
found evidence that the dynamics does not relax onto a quasi-periodic KAM torus, as it has been 
conjectured [niis], but continues to spread chaotically along the chain for arbitrarily long times. 


4 Chaotic diffusion in chains of coupled symplectic maps 


Typically, Hamiltonian systems and in general discrete time dynamical systems can be analyzed 
with the help of discretization techniques, such as Poincare sections. The disadvantage is that the 
resulting map is, in general, not explicitly known and thus one has to resort to the construction of a 
suitable discrete system that exhibits similar statistical properties with its disordered Hamiltonian 
counterpart. Moreover, since the numerical integration of Hamilton’s equations for long particle 
chains is often too costly, one might wonder whether it is possible to use a system of symplectic maps 
to model the dynamics and statistics of multi-degree of freedom Hamiltonian systems, especially in 
the physically interesting case where nonlinearity and disorder coexist [lailllElIlg]. This would 
mean that each particle in such a chain corresponds to a 2-degree of freedom oscillator and is 
modeled by a 2-dimensional symplectic map. 

We decided therefore to describe each particle by the 2-dimensional McMillan map [15] given 
by 

„ , „ „ _ „ , ‘^KXn 

^n-t-1 — T Pn-t-1) Pn +1 — Pn “t” 9,1 2 Xyj. (ff) 

xi + 1 

This map is integrable, possessing the constant of the motion /„ = x^y^ + x'^ + y'^ — 2Kxnyn, 
with In = In+i for all n = 0,1, 2,... and yn = Xn — Pn- If \K\ < !> the origin in the (a;„,p„) 
plane is elliptic, while it becomes a saddle point via a pitchfork bifurcation for K > \ and via 
a period-doubling bifurcation for AT < —1. In Ref. [28], the authors studied a non-integrable 
conservative perturbation of map m and chaotic orbits were investigated from a statistical point 
of view, similar to the approach used here. 

In our study, we shall consider a system of CMMs expressed in the form 


Jj) 

''n+l 


— rrU) 
- 


■Pi+1> 


_ „(j) 

Pn+l ~ Pn 


dV 


dx. 


U) 


j = l,...,TV, 


( 12 ) 


where Xn', pn denote the position and momentum of the j-th map respectively, at the n-th 
step in discrete time (i.e. iteration). Thus, CMMs represent a discrete non-linear system of N 
2-dimensional symplectic McMillan maps with a quadratic nearest neighbor interactions potential 
of the form 

N 

V = ln(xb)2 + 1) + x^^^^ + |(4^+b _ a;^)®), (13) 

under fixed boundary conditions = 0, while e 0 is the coupling parameter that 

renders the system non-integrable. 

In the uncoupled case (i.e. e = 0), if K^^'> > 1, j = 1,... ,TV each map has a saddle point 
at its origin, while for 0 < < 1 the origin is elliptic. Choosing larger than 1, if e > 0 

is small enough, we make sure that the coupling does not alter the stability properties of the 
origin, while the dimensionality of its stable and unstable manifolds can still be adjusted by only 
varying the values of the A'^T) parameters. For these reasons, we have chosen to use throughout 
the paper A'b) values in [1.2,1.6] to emulate disorder as we have additionally checked by numerical 
simulations that any other choice of A'*^^) > 1 ^ produces similar results. Finally, a preliminary 
study for different initial conditions (i.e. initial excitation strengths) revealed that as long as they 
are chosen sufficiently close to the origins of all maps, the obtained results are similar. This is 
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Figure 2: Results for the case of 5 CMMs with e = 0.1. For each row, from left to right, 
we plot the corresponding two largest Lyapunov exponents Ai, A 2 (left panels), the MSD (central 
panels) and qi (right panels) as a function of n. The first row corresponds to = 1.6 and the 
rest of K's equal to 0.6, the second row to the case = 1.6 with all other K's equal to 

0.6, the third panel to the case = 1.6 with the remaining K's equal to 0.6, the 

forth row to the case = 1.6 and = 0.6, and finally the fifth row 

to the case where = 1.6, j = 1,..., 5. In all cases, qi corresponds to ni = and M = 10^. 


supported by the fact that for > 1 all origins are saddle points with stable and unstable 
manifolds and thus a typical initial excitation strength of the j-th map x^'^ = 10“® suffices to 
ensure that the above properties hold. 

Now, if one starts with all maps having stable origins (all < 1) and excites only the central 
map, just as we did in Ref. [3] for single-site excitations of the KG system, no chaotic diffusion is 
observed. We, therefore, turn to the more challenging scenario where the central map has a saddle 
point at its origin and use an ensemble of fVjp initial conditions near that point, assuming that 
all other maps have elliptic points at their origin. Finally, we increase the number of maps with 
saddle point origins and study the dynamics and diffusion of orbits in phase-space by computing 
the corresponding pdfs of sums of their position coordinates Xn as described in Sec. [2]to understand 
if there is any connection between chaotic diffusion in the two models. 


7 























































We first considered N = 5 coupled McMillan maps and studied the dynamics around the origin 
when the central map (i.e. j = 3) is unstable. In the calculations of the pdfs, we used as an 
observable variable rji = of the central map and averaged over iVjp = 10 ^ initial conditions 
randomly chosen within [—10“^°, 10“^°] about a reference orbit started close to the origin. For e 
as small as 0.001 we observed no diffusion since the corresponding 7 exponent was practically 0 
and the pdfs of the sum variable Sm = Sti of the central particle were well approximated 
by q-Gaussians, whose index fluctuates around qi = 2 even after as many as 10® iterations! 

On the other hand, increasing the coupling strength to e = 0.1 allows for the presence of 
stronger chaos and yields signs of diffusive behavior. In Fig. [51 we show for the N = 5 case (with 
M = 10^), in each row, from left to right, the two largest Lyapunov exponents Ai, A 2 (left panels), 
the MSD as a function of the number of iterations n (central panels) and the qi entropic index 
of the central map (i.e. ni = as a function of n (right panels). Note that, as the number 
of unstable maps increases, qi approaches 1 with increasing n and a diffusion process is observed 
whose exponents 7 for large enough n are clearly smaller than 1. In particular, when more than 
two maps are unstable, there is first a transient state of rapid diffusion with MSD oc rC^ and 7 > 2, 
faster than ballistic, over a relatively small number of iterations. Then, as n increases, a second 
epoch occurs with 7 < 1 , corresponding to subdiffusive dynamical processes. 



Figure 3: Characteristic pdfs for 30 disordered CMMs for single site excitations. Panel 
a): Plot of the computed pdf (dashed) for the observable ryis = (that represents the 

15 central maps) for the time interval [0,10®] with qi^ = 1.38 ± 0.05, taken from fitting with a 
g-Gaussian distribution © in solid thick. Panel b): Similar plot of the computed pdf (dashed 
curve) for the observable 730 = Si/=i maps) computed for the time interval [ 0 , 10 ®] with 

gao = 1.04 ± 0.05, taken from fitting with a g-Gaussian distribution © in solid thick. In both 
panels, e = 0.28, N = 30, M = 5 • 10® and K values were chosen randomly in [1.2,1.6] to 
implement disorder. Both cases were taken from the second panel of Fig. [S] Note that the vertical 
axes are in logarithmic scale, while the dotted curve is the Gaussian pdf with g = 1 to guide the 
eye. 

The next step is to set TV = 30 and increase the coupling further setting e = 0.28 with single site 
excitations perturbing only the central map, and all values chosen randomly in the interval 
[1.2,1.6] to emulate disorder as in the KG system. In this case, we obtain much more interesting 
statistics. As can be seen in Fig. O strong chaos is displayed when all maps are taken into account 
(panel b)) with gag = 1.04 ± 0.05 (that corresponds to 7730 = while the 15 central 

maps of the chain exhibit weak chaos (panel a)) with gi 5 = 1.38 ± 0.05 (where 7715 = 

These findings are important as they demonstrate a similar situation to the physical problem of 
disordered KG coupled oscillators studied in Ref. [3], where the central particles perform weakly 
chaotic motion, while the overall system displays strongly chaotic dynamics depicted by their 
corresponding entropic indices. 

Based on these first results, our next step was to investigate CMMs of increasing e to relate 
complex statistics with diffusive dynamics. We did that for 30, 50 and 100 disordered CMMs using 
single site excitations. In all cases studied, we witnessed the same scenario: In order to observe 
(sub)diffusion in the disordered CMM, the coupling strength had to be relatively large (e > 0.05) 











Figure 4: Results for the case of 50 disordered CMMs and single site excitations with 

e = 0.1. Panel a): Plot of the MSD as a function of n for = 2 • 10^, all maps locally 
unstable by using K values randomly chosen in [1.2,1.6] with only the central map initially excited. 
Panel b) shows the evolution of the q entropic index as a function of n for the central map {qi 
with rji = the 25 central maps (925 with 7725 = ^0 maps (550 for 

?750 = X]j=i considered. The computed pdf (dashed curve), the fitted g-Gaussian (solid thick 
curve) of Eq. 0 and the Gaussian pdf (dotted curve) are shown at n = 10® iterations: Panel c) 
for the central map, d) for the 25 central maps and e) for all 50 maps. In all panels for the entropic 
index g, we have used M = 5 • 10®. 

with all maps being locally unstable at the origin with all randomly distributed in [1.2,1.6]. 
We thus concluded that, for small enough coupling, orbits do not diffuse in phase-space but rather 
remain wandering around the unstable origins of each GMM. This behavior was clearly depicted 
by the entropic indices as the central particle performed a weakly chaotic motion with qi values 
substantially higher than unity whereas the rest of the chain behaved in an essentially strongly 
chaotic manner approaching g ft! 1 in the long run. 

We then turned to bigger values of e, focusing on 30 and 50 CMMs as the case of 100 GMM 
performed similarly. In particular, for 50 CMMs and e = O.I, exciting only the central map and 
taking all maps to have a saddle point at their origin, we found that after an initial number of about 
10^ iterations (during which the orbits diffuse very fast with MSD oc the chain settles down 

to a clearly subdiffusive regime with MSD oc n® ® as evidenced in the hrst panel of Fig. [H This 
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behavior, however, is associated with weak chaos, since the entropic indices of panels b) to e) of the 
same figure, for the central map {qi for 771 = the 25 central maps (525 for 7725 = Sjlia 

and the whole chain (950 for 7750 = ) as high as 2 for large nl 

Using then coupling strengths as high as 0.28 for 30 and 0.29 for 50 CMMs, presented in Fig. 0 
we did observe fully chaotic behavior for all different parts of the chain accompanied by subdiffusive 
motion. This is illustrated in the first two panels of this figure for 30 CMMs and the last two for 50 
CMMs. Clearly, in the long run, the subdiffusive motion is characterized by an exponent 7 ~ 0.68 
and by all entropic indices fluctuating in the interval [1,1.2], signaling that the dynamics of the 
chain is close to fully developed chaos. These findings provide evidence that the CMM system in 
the long run, for sufficiently high dimensionality and coupling strength,sup does show evidence of 
strong chaos and subdiffusion, and thus exhibits dynamics quite similar to what was observed in 
the KG continuous flow model studied in Ref. [3]. 





Figure 5: Subdiffusive chaotic behavior of large CMM disordered systems {N — 30 and 

50). Panel a): Plot of the MSD as a function of n ioi N = 30 CMMs, e = 0.28, iVjj, = 4 • 10^, 
all maps locally unstable and K values randomly chosen in [1.2,1.6] with only the central one 
initially excited. Panel b) shows the evolution of the q entropic index as a function of n for the 
central map (gi for 771 = the 15 central ones (<715 for 7715 = maps ((730 

for 7730 = when considering TV = 30 CMMs. Panels c) and d) are similar for N = 50, 

e = 0.29 and = 2 • 10^. Note that all horizontal and vertical axes of the MSD plots are 
logarithmic. In all panels for the entropic index g, we have used M = 5 • 10®. 


5 Conclusions 

In this paper we investigated the connection between regimes of “weak” and “strong” chaos and 
diffusive dynamics in a multi-dimensional system of coupled symplectic McMillan maps, with the 
aim of comparing our results with analogous studies by S. Flach and coworkers on a Klein-Gordon 
Hamiltonian in the presence of disorder and nonlinearity. Our main idea was to explore whether 
coupled symplectic maps can be used to study wavepacket subdiffusive spreading in disordered 
media, thus avoiding the tedious integration of large systems of ordinary differential equations that 
describe Hamiltonian particle chains such as the KG model. 
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Our statistical approach was based on the computation of sums of position coordinates, in 
the spirit of the central limit theorem and nonextensive statistical mechanics, approximating the 
probability distribution functions by g-Gaussians, whose index g > 1 is connected with weak chaos, 
and g = 1 with strong chaos. 

Let us summarize our findings in reference to the questions we raised in Sec. [H First, we 
demonstrated that our CMM analogue of a disordered nonlinear system of Hamiltonian oscillators 
displays chaotic diffusion connected with two main epochs: A relatively short one of super-diffusion 
followed by a second subdiffusive stage with exponents quite smaller than 1. This is in contrast 
to diffusive motion in a disordered KG chain where only one epoch is observed [ISI3- Secondly, 
we studied GMMs with increasingly higher dimensionality to explore the relation between complex 
statistics and subdiffusion and verified in all cases that to observe subdiffusive motion the coupling 
strength has to be moderately large and independent of dimensionality. Thus, in comparison with 
our findings for the KG system in Ref. [5], the dynamics and statistics of diffusive motion in the 
CMM model yield analogous results. 

Finally, we showed that subdiffusion is accompanied by weak and strong chaos (with g > 1 and 
g = 1-Gaussian pdfs respectively) regardless of the dimensionality of the system, but depending 
on the coupling strength, in analogy to a disordered KG chain at different energy ranges. The 
key point behind this remark is that in the KG system nonlinearity is measured by the amplitude 
of the initial excitation that determines the total energy and diffusive/statistical properties of the 
dynamics, whereas in the coupled map model chaotic diffusion is measured by the magnitude of 
the coupling strength which is proportional to how much the GMMs deviate from the integrable 
e = 0 limit (see Eq. 0). Finally, in the coupled map model disorder does not seem to play an 
equally important role as in the KG system. 

Thus, even though there are substantial differences in the local dynamics of the two systems 
considered here, our results suggest that coupled symplectic maps may be used to model subd¬ 
iffusion and weakly chaotic behavior in certain classes of multi-dimensional nonlinear oscillators. 
For such systems, employing discrete iterative systems can circumvent costly numerical algorithms 
needed to integrate the corresponding continuous time equations of motion and identify maps 
through complicated discretization techniques. 
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